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We simulate magnetic resonance force microscopy mea- 
surements of an entangled spin state. One of the entan- 
gled spins drives the resonant cantilever vibrations, while 
the other remote spin does not interact directly with the 
quasiclassical cantilever. The Schrodinger cat state of 
the cantilever reveals two possible outcomes of the mea- 
surement for both entangled spins. 

Magnetic resonance force microscopy (MRFM) pro- 
posed a decade ago Q is now approaching its ultimate 
goal: single spin detection The following question 
arises: To what extent can the MRFM be used for quan- 
tum measurement of spin states? The particular problem 
considered in this paper is the MRFM measurement of 
entangled spin states using a cyclic adiabatic inversion 
of the spin, which drives the resonant vibrations of the 
cantilever. We discuss the possibility of determining the 
state of a remote spin which is entangled with the spin 
interacting with the MRFM measurement apparatus. 
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FIG. 1. Diagram of an MRFM measurement of an entan- 
gled state of two spins in a chain of spins. Si is a single mea- 
sured spin which changes its direction under the action of the 
alternating magnetic field B\ . §2 is the remote spin entangled 
with the spin Si. m is the magnetic moment of a ferromag- 
netic particle which is attached to the cantilever tip. The 
magnetic force on the cantilever tip (attractive or repulsive) 
depends on the direction of the spin Si . Bo is the permanent 
magnetic field which points in the positive z-direction. 

As an example, we consider a typical entangled state 



for two spatially separated spins, 

(l/^/2)(| TT) + Ui)). 



(1) 



According to the conventional point of view, by measur- 
ing the left spin in the state | t) (or | i)) we automatically 
collapse a remote right spin into the same state | t) (or 
I D). In the process of MRFM measurement (see Fig. 1) 
the direction of the measured spin. Si, changes periodi- 
cally with the period of the cantilever vibrations. Thus, 
it is not clear what the direction of the entangled remote 
spin, 5*2, will be after the MRFM measurement. 

To study this problem, we simulated the quantum dy- 
namics of this spin-cantilever system assuming that the 
measuring spin is initially entangled with the remote 
spin. The remote spin is not subjected to the action 
of the MRFM apparatus. 

The dimensionless quantum Hamiltonian of the spin- 
cantilever system in the rotating reference frame is [p|. 



(2) 



Here Pz and z are the dimensionless momentum and co- 
ordinate of the cantilever tip; Si is the "first" measured 
spin; e = eir) is the dimensionless amplitude of the radio- 
frequency (rf) field (where r = LOct is the dimensionless 
time and Uc is the cantilever frequency); rj is the dimen- 
sionless constant of interaction between the cantilever 
and the spin, which is proportional to the magnetic field 
gradient produced by the ferromagnetic particle on the 
cantilever tip. The phase of the rf field is taken in the 
form {cot + v(i))j where to is chosen equal to the Larmor 
frequency of the spin: lo = ujl. The time derivative, (p, 
changes periodically with the frequency of the cantilever 
vibrations, lUc. In our notation, the dimensionless fre- 
quency of the cantilever vibrations is one unit. Thus, 
the dimensionless period of the cantilever vibrations is 
2tt. The periodic oscillation of (f provides a cyclic adia- 
batic inversion of the spin [|j , which drives the resonant 
vibrations of the cantilever. 

The dimensionless wave function of the whole system, 
including the second entangled spin, can be written in 
the ^-representation as, 

^{z,Szi,Sz2,t)^u^^{z,t)\ tt) +MTi(^.^)l Ti)+ (3) 
uu{z,t)\ i^)+Uii{z,T)\ ii). 
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Substituting (3) into the Schrodinger equation, we derive 
four coupled equations for the functions r), 



(4) 



This system of equations sphts into two independent sets 
of equations. 




FIG. 2. The probability distribution, P(z), at r = 216. 

The initial condition is assumed to be a product of 
the coherent quasiclassical wave function of the cantilever 
and the entangled state of the two spins, 

*(z,&i,^.2,0) - (1/a/2K(z)(| TT> + I ii)), (5) 



Ua{z) = TT-i/^ exp[-(z - V2af/2], 

where zq = \/2a is the average initial coordinate of the 
cantilever. For numerical simulations we used the follow- 
ing values of parameters, 

a = -IOV2, 77 = 0.3, 

and the following time dependences for e(T) and (p{t), 

e = 20t, t < 20; e = 400, r > 20, (6) 

93 = -600 + 30r, r < 20; ip = 1000sin(T - 20), r > 20. 

The chosen value of a corresponds to the quasiclassical 
state of the cantilever with the number of excitations 



n = |ap = 200. The chosen value of rj corresponds to 
the existing MRFM experiments with electron spins 
The time dependences (6) provide the conditions for the 
cyclic adiabatic inversion of the spin. Fig. 2 shows the 
typical probability distribution of the cantilever position. 



P(z) = |unl' + l«TiP + l"iTp + l"ul' 



(7) 



One can see that the probability distribution, P{z), de- 
scribes a Schrodinger cat state of the cantilever with two 
approximately equal peaks. When these two peaks are 
clearly separated, the total wave function can be repre- 
sented as a sum of two terms corresponding to the "left" 
and the "right" peaks in the probability distribution, 

^iz,Szi,Sz2,T) = *q(2:,S'^i,S'22,t) +-^b{z,Szi,Sz2,T). 

(8) 

Our numerical analysis shows that each term, and 
can be approximately decomposed into a direct product 
of the cantilever and spin wave functions, 

^'a = Ua{z,T)Xa{Szl,T)\ t)2, ^'h = Ui,{z , T)xb{S zl, t)\ 1)2. 

(9) 

This decomposition is possible because the complex func- 
tion U||(z,r) is proportional to U||(z,r), and the com- 
plex function M||(z,t) is proportional to uii(z,t). Such 
proportionality can be seen in Fig. ^, where we plot the 
corresponding wave functions at the same time as in Fig. 
2, with suitable numerical coefhcients. 
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FIG. 3. Upper boxes: the real part of wave functions; 
lower boxes: the imaginary part of wave functions, at r = 216. 
a)iie(it||): solid line; i?e(— 5u||): filled circles, b) Re{u^i): 
solid line; i?e(5M||): filled circles, c) Im{u]]): solid line; cir- 
cles /m(—5Mj|)filled circles, d) 7m(u||): solid line; /m(5«|j.): 
filled circles. 

The spin wave function, Xa{Szi,T), describes the dy- 
namics of the first spin with its average, (Xa|'5'|Xa), point- 
ing approximately in the direction of the effective mag- 
netic field, Bes, in the rotating frame. 



Bcs = (e, 0,-c^). 



(10) 
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(We neglect here the nonhnear term 2riz whose contribu- 
tion to the effective field is small.) 
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FIG. 4. The positions, Zm,ax(j), of two peaks of the 
Schrodinger cat state as a function of time. 

The spin fimction, ^^(S'^Ijt), describes the dynamics 
of the first spin with its average pointing in the direction 
opposite to the direction of -Boff- As the amplitude of 
the cantilever vibrations increases, the phase difference 
between the oscillations of the two peaks, \ua[z^ r)p and 
|ufc(z,T)p, approaches tt. (See Fig. 4.) (To reach the 
phase difference of tt, a long time of numerical simulations 
is required. That is why we restricted the simulation time 
for the results presented in Fig. 4.) 




FIG. 5. Two outcomes of the MRFM measurement of the 
state of two entangled spins, (a) The measured spin, Si, 
points along the direction of the effective magnetic field, and 
the remote spin, 5*2, points "up" (in the positive z-direction) . 
(b) The measured spin. Si, points in the direction opposite to 
the effective magnetic field, and the remote spin, S2, points 
"down" (in the negative 2-direction) . 



In realistic experimental conditions, the Schrodinger 
cat state quickly collapses due to the interaction with 
the environment In this case, the two peaks of the 
probability distribution describe two possible trajectories 
of the spin-cantilever system. In one of these trajectories 
the first (measured) spin is pointed along the direction of 
the effective magnetic field while the second (remote) spin 
is pointed "up" (in the positive z-direction); the other 
trajectory corresponds to the opposite situation in which 
the orientation of both spins is reversed — the first (mea- 
sured) spin is antiparallel to the effective magnetic field, 
and the second (remote) spin is pointed "down" (in the 
negative z-direction) . The phase difference between the 
corresponding oscillations of the cantilever approaches tt 
with increasing cantilever vibration amplitude. 

In summary, we have studied the outcome of the 
MRFM measurement of the entangled spin state, 
(I/\/2)(| tt) + I ii))- Our numerical simulations reveal 
two possible outcomes shown schematically in Fig. 5: 
(a) The first (measured) spin points along the effective 
magnetic field in the rotating frame, and the second (re- 
mote) spin points in the positive z-direction. (b) The 
first (measured) spin points opposite to the direction of 
the effective magnetic field in the rotating frame, and the 
second (remote) spin points in the negative z-direction. 
Thus, the collapse of the measured spin along (or oppo- 
site) the direction of the rotating effective magnetic field 
leads to the collapse of the remote spin in the positive (or 
negative) z-direction. These two outcomes correspond to 
two phases of the cantilever vibrations which differ by tt. 
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